Effect of light cr-meson Production in pp — > 37r° at rest 
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The 7r°7r° mass spectra and angular distributions around A'A'-threshold and at 1.5 GeV 
in pp(at rest) — > 3tt° in the Crystal Barrel experiment are reanalyzed by applying the new 
method, which is consistent with unitarity of S-matrix and expressed directly by resonance 
parameters. The effects of light cr-meson production are clearly seen to improve the fit with 
a, in comparing with the fit without a. 

§1. Introduction 

The light iso-singlet scalar a meson plays an important role in the mechanism of 
spontaneous breaking of chiral symmetry, and to confirm its real existence is one of 
the most important topics in hadron physics. Recently in various 7T7r-production ex- 
periments 1 " 2 ) a broad peak is observed in mass spectra below 1 GeV. Conventionally 
this peak was regarded as a mere non-resonant background, basing on the "univer- 
sality argument," 3 ) since no a was seen in tttt scatering at that time. However, at 
present the 7T7r-scattering phase shift 5g =0 is reanalyzed by many authors 4 )> 5 )' and 
the existence of a meson is strongly suggested. A reason of missing a in the conven- 
tional analysis is pointed out to be due to overlooking the cancellation mechanism, ^ 
which is guaranteed by chiral symmetry, between the effects of a and those of repul- 
sive 7T7r-interaction. Moreover, the conventional treatment, based on the universality 
argument, of the low mass broad peak was shown not to be correct, and a new ef- 
fective method, variant mass and width(VMW), is proposed to analyze resonance 
productions. In this method the production amplitude J- is directly represented by 
the sum of Breit-Wigner amplitudes with production couplings and phase factors 
(, including initial strong phases) of relevant resonances. The consistency of this 
method with the unitarity is seen from the following field theoretical viewpoint ^ . 

Presently, after knowing the quark physics, the strong interaction £ str among 
hadrons (in our example, mesons) is regarded as a residual interaction of QCD 
among color singlet (/(/-bound states, the "bare states," denoted as ir = ff , a, fo and 
/2- In switching off C s t r , the bare states appear as stable particles with zero widths. 
In switching on £ s tr, they change into the physical states with finite widths. The 
unitarity of S-matrix is guaranteed automatically by the hermiticity of £ s tr- The 
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VMW method is obtained directly as the representation of production amplitude by 
physical state bases, with a diagonal mass and width. 

The VMW method is already applied 1 -* to the analyses of pp-central collision 
pp — > ppTT°TT° and J/ip — > uuinr decay, leading to a strong evidence of existence of the 
light a meson. In this paper we apply this method to analysis of the high statistics 
data on the process pp — > 3ir° at rest obtained in the Crystal Barrel experiment. 2 ) 

§2. Amplitude describing pp — > 3ir° at rest by VMW method 

2.1. £ str for pp — ► 3-7T at rest 

We apply the iso-bar model, describing the process in following two steps: In the 
first step the pp annihilates into the resonance f(fo and $2) and 7r°, and the / decays 
into 27T° in the second step. Since both p and p are at rest, the relative momentum 
Vii = Ppfi—Ppfj, {and the relative angular momentum Lpp) between p and p is 0. Charge 
conjugation parity Pc of /7r°-system is +1. Thus the three types of (p,p) bi-linear 

forms with Pc=+1 are possible: pi^5P,pi^5^ fJ ,p and pp. The second type reduces to 

<— ► 

the first type by using the equation —d^pi^p) = 2m p pi^^^p + pi^io^ d u p = 
2m p pi^^ lJL p (, derived from Dirac equation and the above "rest condition"), and 
the third type is forbidden by parity. Thus only the first type remains.*) The most 
simple form of describing the 1st step pp — ► fir and the 2nd step / — > 2ir° is 
given, respectively, by 

£str = J2_ (^/o^75P/07r° + £/ 2 pi75P/2^<9/A7T ), 

F0J2 

£lr = E (9f fW + 5/ 2 /2^(vr tt)). (2-1) 
/0J2 

2.2. Amplitude by VMW method 

First denoting the three 7r° as 7Ti,7T2 and 7r3 with momenta p±,P2 and j>3, re- 
spectively, we consider the 7Ti and iT2 forming the resonance / with squared mass 
s±2 (where Sij = —(pi + Pj) 2 ) and with momentum |p| in z-direction. The tt\ has 
momentum |q| and polar angle 9 in the / rest frame. In the lowest order in bare 
state representation 6 ^ the amplitude is **) 

pSp( ST iMlo , V Zh9f 2 N( Sl 2, cosM , 



2im p f^(^ + ^ -), 

where 



- S12 V m fo ~ s i2 

/() /0 /2 /2 



2im p f SpSp = p(0, Sp)i75p(0, s p ) (s p and s p being spin of p and p) 



*•* This implies that the initial pp is in the 1 5'o state. However, in the original analysis 2 ', the 
phenomenological parameters related with the Pi and P2 states, which are considered not to 
contribute since of the above "rest condition," are included. 

**) N(si2,cos 12 ) is obtained by the calculation of N = —pz^Pz v P^ v -\ K {pi - P2)\(pi - P2) K , 
where we use the tensor projection operator V^Xk with mass squared S12 instead of rn 2 f 2 : P^ V -\ K = 
\{&n\5 VK + <5 mk (5„a) — \8^ V 8\ K , where <5 M! , — 8^ v + F ^ ; P M = pi M +p2 M . 
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ai(Si2 , cos ,, ) = _(£21^ + 16 !!!|pV 
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Owing to the effect of final (and initial) state interaction, the "full order" of the 
amplitude in physical state representation is given by 



A SpSp (s 12 ,cos8 12 ) = 2im p f" s P V — 

\r ml - si2 - i-i/s 



/() no fo - S12 ~ iyfs~^r fo {s 12 ) 

+ y r /2 e^Ar( gl2 ,cosg 12 ) \ ^ 

The symmetric amplitude, satisfying the statistics property of 3tt system, T s 8f is 
obtained simply by its cyclic sum as 

F Sp s p = A SpSp (s 12 ,cos9 12 ) + A SpSp (s 23 ,cos6 23 ) + A SpSp (s 31 , cos6» 3 i ) . (2-4) 

Cross section is given by 

r2m p -m n Hf^ | p | | q | rl 

da ~ / dy/s^^ — — / dcosO \F\ 2 , 

Jim-* vr 87rm p stt^/s^ J-i 

where 

^^(l/4)^|^ pSp f = (l/2)|^ + _| 2 , 

s\2 + s 23 + S13 = 4m 2 + 3m 2 , 

s 23 = -4m p (|p||q|/^/sl2")cos6l + m 2 + 2m p E 3 , 

E 3 = ym^+p2 = ( 4m 2 _ + m 2^ 4j7lp _ (2-5) 

§3. Results of Analysis 

We analyze the experimental data, the 7r°7r° mass spectra and angular distri- 
butions around JT^-threshold and at 1.5 GeV, which are published in the paper 2 ) 
by Crystal Barrel collaboration, using our formulas Eq.(2-5). We take into con- 
sideration as the physical particles /o = a, /o(980), /o(1370), /o(1500), and f 2 = 
/ 2 (1275), / 2 (1565). 

We take into account the tttt and KK couplings of the relevant resonances, 
where we consider the effects of all the inelastic channels are represented by the KK 
coupling. 

The result of the fit is shown in Fig. 1. 

The mass and width of a obtained are m a = 540 i^MeV and T a = 385 ^soMeV 
(error corresponding to the 5<r-deviation) . The reduced \ 2 is given by 746.9/(281 — 
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Fig. 1. Result of the fit by the amplitude represented in terms of the resonances with tttv and KK 
couplings: (a) 7r°7r° mass distribution, and cos# distributions (b) around if A'-threshold and (c) 
at 1.5 GeV. Theoretical curves of the cos# distributions around Tfif-threshold and at 1.5 GeV, 
are given at the energies, 995.3 MeV and 1500 MeV, respectively. The total \ 2 is Xtotai = 746.8, 
as a sum of the respective contributions from (a), (b) and (c), 406.5 + 148.1 + 192.2. The reduced 
X 2 is Xtotai/(^data ~ -/V param ) = 746.9/(281 — 30) = 2.98. Dashed curves represent the spectra 
obtained by setting r CT = in the fit. 

30) = 2.98. Respective contributions to the x 2 from the mass and angular distri- 
butions(a.d.) around KK and at 1.5GeV are 406.5(mass), 148.1(a.d. around KK) 
and 192.2(a.d. at 1.5GeV).*) **•* The properties of all the relevant resonances used 
in the fit are given in Table 1. 

We have also tried the fit without considering the KK couplings of the reso- 
nances. The result of our fit with the KK couplings is almost the same as the one 
without the KK couplings; and in order to determine the KK couplings and the 
other ones, it is necessary to analyze the data on the corresponding channels directly. 

In order to see the effect of cr-meson production in our fit, in Fig. 1 the spectra 
given by setting r a = are also given by dashed lines. Effect of cr-production is seen 
to be crucially important in reproducing the structure of mass spectra below 1 GeV. 

We have also tried the fit without introducing the a Breit-Wigner amplitude. 

*-* The obtained x 2 value may be compared with the one by the original fit 2 ': We have estimated 
their corresponding x 2 by reading the deviations of their fit from the experimental points and 
obtained the numerical values as x 2 /{Ndata ~ N para m) = 835.4/(282 — 25) = 3.25. Respective 
contributions to the \ 2 are 432.8(mass), 113.4(a.d. around KK) and 289.3(a.d. at 1.5GeV). Almost 
the similar but slightly improved fit is obtained in the present VMW method. However, note that 
their parameters were determined by analyzing the original two-dimensional Dalitz plot directly. 
They reported that \ 2 in this case is /N F = 2028/(1338 - 34) = 1.6. 

**' Number of data points in the figures given in ref.2) is 282=82(mass) + 100(a.d. around 
KK)+W0(a,.d. at 1.5GeV). In our analysis the first one point close to the threshold of mass spectra 
is removed since this point is at 279MeV, where -k 7T channel is open but ir + ir~ channel is almost 
closed, and accordingly a special treatment of -kit widths of resonances is required. 
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Table I. Values of the resonance parameters obtained by the \ 2 fit with both tttt and KK 
couplings^TTTr, §kk)- The errors correspond to 5 a deviation. The mark * represent that 
the corresponding values are not able to be determined by the fit: The gKK(a.nd corre- 
spondingly Fkk) of /2(1270) is insensitive to the fit. The mass of /2(1270) falls in its con- 
strained lower limit, and the tttv(KK) coupling of /o(980) does in its upper(lower) limit. In 
our analysis energy- dependent widths of the resonances are given by the sum of the tt-k and 
KK widths, r(s) = f„(s) + f K ^(s), where /^(s) = 3i gf s 1 (i = nn,KK) for scalars and 

4 „2 |p 1 5 _ 

-T;(s) = {i = tttt, KK) for tensors. The values of _T in the table are the widths at 

the respective peak positions of the resonances: r = r(s = mj es ), except for the F KR of 
/o(980) (, of which value is calculated as F KR = i f^ ds r KR (s) (m ^ a ^ s) 2 1 +artot(s) 2 \ N = 

f , ds t — 5 — - — rrp)- The relative ratios of the absolute magnitude of angular dis- 

tributions compared with the mass spectra are, respectively, 0.0695 and 0.1606. The cen- 
tral values of production coupling r and of production phase 0(deg.) are given, respec- 
tively, by (r -,r /o (g 80 ),r /o( i3 7 o),r /o (i5oo), ^(1270)^/2(1565)) = (370,108,1720,676,1721,4724) 
and (^(ggo), (9/ (i370) 5^/0(1500)1^/2(1270)1^/2(1565)) = (258, 164, 111,276,71). 





mass(MeV) 


^(GeV) 


r^(MeV) 


tex(GeV) 


r KR (MeV) 





+36 

540-29 


3.47±S:1I 


385±1^ 






/o(980) 
/o(1370) 
/o(1500) 


964±w 

-1-19 

1368II6 
1518±i? 


2.25 
4.10 + -°f 

i.9oi8:2S 


100 

240± 3 * 
46±lf 


2 

* 

2.65+- 10 , 1 


* 

70± 63 


/a (1270) 
/a (1565) 


1220 
1552±2i 


8.6oio:6° 

7.88 + - 3 -f 


142±IS 

gst 12 , 1 


16 

14.7 + - 3 f 


38 

ioo± 5 t 



In this case the broad peak structure below 1 GeV in the mass spectra are only 
roughly reproduced by the combinatorial background coming from the higher mass 
2-7T resonances due to the statistics property of 3-7T system. The corresponding x 2 
is 3112/(281-26)=12.2, which is much worse than our best fit with a meson. This 
seems to give a strong evidence for u-existence. 

§4. Comparison with other analyses 

Several extensive analyses (including the original one 2 )) on the relevant exper- 
imental data by Crystal Barrel collaboration have been thus far done. However, 
all the analyses seem to be done, more or less, under the influence of "universality 
argument 3 )." 

According to this argument, all the irn production amplitudes T must be parame- 
trized through the "universal" tttt scattering amplitude T as 

F = a{s)T, (4-1) 

with slowly varying real function a(s) which corresponds to the nir production cou- 
plings of the production channel and is process-dependent. This equation is believed 
to be based on the unitarity or final state interaction (FSI) theorem. 

However, it has been pointed out that, in order to apply the FSI theorem 
presently, after knowing quark physics, we must make a special attention 6 ) on the 
bases of representation of T and T. Before knowing quark physics, the S'-matrix 
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of strong interaction S s t T was represented in terms of only stable particles such as 
7r and N. However, now S str should be described by the interaction Hamiltonian 
Ti i ,str among the bare states, the color neutral stable bound states of quarks and/or 
antiquaries. The observed irir resonances, such as a or /q(980), and all the resonant 
states must be treated equally to the stable states, such as one pion or two pion 
states, as complete set describing the hadron world. They have mutually indepen- 
dent 7T7r( and/or KK) production couplings, in principle. Accordingly, Eq.(4-1) in 
its original form has proved not to be correct. 

In the pioneering work 7 ) by Aker et al. of Crystal Barrel(CB) collaboration, the 
I = S wave production amplitude Fs{si2) (in our notation) is taken as 

^12) = a —L_e^° (s ^ sin 4 = °(*i2) = aT (4-2) 

with the constant tttt coupling a, basing on the original "universality argument" 
Eq.(4T). The symmetric production amplitude is obtained by taking the cyclic sum 
as Eq.(2-4), while the scattering amplitude T is picked up from the reference by 
Au, Morgan, Pennington(AMP) in 3), where the conventional JC matrix analysis was 
done for the CERN-Munich tttt scattering phase shift 5g =0 . 8 ) 

In a series of works by Anisovich, Sarantsev, Bugg and Zou et al. 9 ) " n ) being 
done in the line of this thought 7 ) , the relation of "the universality argument" T = 
a(s)T with real a(s) function is argued to be not correct, since of the possible 
effect of the strong phases(basing on the detailed consideration of their origins, for 
example, triangle singularities). They used N/D formalism and mentioned that 
in production processes the complex N function N'(s), being independent of -kit 
scattering N function N(s), is necessary. A great variety of modification is allowed 
for the parametrization of N'(s), and they try to fit the spectra of pp annihilation 
with three tentative forms of N'(s) with complex couplings. Here, they introduce 
the scalar (/ (1365) and / (1520)) and tensor(/ 2 (1270) and / 2 (1560)) Breit-Wigner 
amplitudes above 1.1 GeV with complex production couplings. In our interpretation, 
their analysis is, so to speak, two-fold: Below 1.1 GeV it was done without taking 
into account of the freedom of production couplings of resonances by following the 
viewpoint of "the universality argument," while above 1.1 GeV their method is 
equivalent to the VMW method, where this freedom is explicitly introduced. 

In all the above analyses the pole positions of T matrix in low energy region 
below ~ 1 GeV are determined only through the analyses of CERM-Munich 5^ =0 8 ) by 
K. matrix method, although there are many varieties of parametrization methods of 
/C. This is also the case in the original analysis 2 ) by Amsler et al. of CB collaboration. 
In this analysis the tttt scattering and production amplitudes are given, respectively, 
in the K, matrix representation as: T = /C/(l — ipK.), T = V/(l — ipK.), where 
ZC-matrix and "P-matrix are taken in pole-dominative form; 

Here the summation a is taken for the ZC-matrix states, which are related 6 ) to 
the physical states corresponding to the poles of T (or J 7 ). In the case with no 
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production phases, 9 a = 0, the T and T have the same phase, coming from the 
common factor 1/(1 — iptC), and the final state interaction theorem is satisfied. 
Actually in their analysis this phase was set to be the experimental scattering phase 
shift (5g =0 . However, in the above /C-matrix parametrization, when s is close to m 2 , 
K. diverges and the phase must take the value 90°(+nx 180°). This gives a very 
strong constraint for the value of m a . The experimental e)g =0 passes through 90° at 
about \fs ~ 900 MeV, and so the m a becomes m a ~ 900 MeV , which is much larger 
than m (T (~600MeV). Thus, no existence of light a is implicitly assumed from the 
beginning. This situation is common in all the analyses 7 )' 9 ) " n )' 2 ) thus far made. 

In our method, the (5 S =0 is analyzed by introducing the repulsive background 
phase shift <5bg> which is required from chiral symmetry. 12 )' 13 )' ^ The scattering 
S'-matrix and correspondingly the /C-matrix are parametrized by 

S = S *<"S BG , JC= ^j-^G (44) 

1 - p 2 /C Rcs /C BG 

The /C Res in denominator of /C removes the poles of /C Res = J2 a ga/( m a ~ s ) i n the 
numerator in the total /C-matrix and we can take the light m a ~ 600MeV.*) 
Correspondingly, the T is represented by **^ 

•pRes 

which satisfies the final state interaction theorem. This T is able to be 6 ^ rewritten 
into the form, applied in VMW method, in the physical state representation. In our 
approach, whether the light a meson exists or not is determined directly from the 
experimental data themselves, as was done in §3. 



§5. Conclusion 



Through the results of analyses given above we may conclude that the effects 
of production of light a meson are clearly shown. The numerical values of mass 

and width of a are obtained as m a = 540 -29MeV, F a = 385 ^loMeV, which 
are consistent with those obtained in our phase shift analysis 5 ) ((m<j, r , cr )=(535 ~ 
675, 385 ± 70)MeV). However, the effect of a in this process is, in principle, not able 
to be discriminated from those of higher mass resonances, such as /o(1370), due to 
the effects coming from statistics property of 3tt° system. In order to avoid this, it 
is desirable to analyze also the process, pn — > 7r°7r°7r~, through the similar method. 

Finally it should be noted that the experimental data of the spectra applied in 
this paper were obtained by reading out the corresponding figures of ref. 2) and 
incomplete. The excellent reproduction of the data encourages us to study more in 
details. It is desirable to reanalyze directly the experimental data of the Dalitz plot. 

*^ On the other hand, the background matrix, cbg, in the conventional /C-matrix cannot describe 
the global phase motion corresponding to Sbg in our method, and the small value of m a is not 
permissible. 

**' Here we neglect the possible effect of non-resonant 37r° production. 
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